Mon. Not. R. Astron. Soc. 0Q0,[Tlfl9l() Printed 4 May 2009 (MN KT^ style file v2.2) 



The globular clusters-stellar haloes connection in early type galaxies 
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ABSTRACT 

This paper explores if, and to what an extent, the stellar populations of early type galaxies 
can be traced through the colour distribution of their globular cluster systems. The analysis, 
based on a galaxy sample from the Virgo ACS data, is an extension of a previous approach 
that has been successful in the cases of the giant ellipticals NGC 1399 and NGC 4486, and 
assumes that the two dominant GC populations form along diffuse stellar populations sharing 
the cluster chemical abundances and spatial distributions. The results show that a) Integrated 
galaxy colours can be matched to within the photometric uncertainties and are consistent with 
a narrow range of ages; b) The inferred mass to luminosity ratios and stellar masses are within 
the range of values available in the literature; c) Most globular cluster systems occupy a thick 
plane in the volume space defined by the cluster formation efficiency, total stellar mass and 
projected surface mass density. The formation efficiency parameter of the red clusters shows a 
dependency with projected stellar mass density that is absent for the blue globulars. In turn, the 
brightest galaxies appear clearly detached from that plane as a possible consequence of major 
past mergers; d) The stellar mass-metallicity relation is relatively shallow but shows a slope 
change at A/* « IQ^^Mq. Galaxies with smaller stellar masses show predominantly unimodal 
globular cluster colour distributions. This result may indicate that less massive galaxies are 
not able to retain chemically enriched intestellar matter. 
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1 INTRODUCTION 

' Pioneering ideas about the role of globular clusters (GC) 
as tracers of the formation history of galaxies come from 

lEggen. Lvnden Bell & S andag d ( 1 96211 or. with a different perspec- 
tive, from lSearle & Zina,(1978ll . The ev olution of this subject dur - 
ing recent years has been reviewed in ISrodie & StradeJ j2006l) . 
A key issue in that context is the connection between GC and 
the "diffuse" stellar populations of the galaxies they are associ- 
ated with. Some reasons allow to infer that such a link might not 
be simple. However, several papers point out that massive early 
type galaxies seem to ha ve highly synchronised ages. For example, 

iLa Barbera et alj ^20081) reach that conclusion in their discussion 
of the properties of the Fundamental Plane for early type galax- 
ies also indicating that the age variation per mass decade seems 
smaller than a few percent. This suggests that the eventual connec- 
tion between field stars and GC may still be "read" since the effects 
of subsequent (younger) stellar populations in these galaxies may 
have not been important. 

The theoretical side of the problem has received a num- 
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ber of contributi ons in an effort to b etter define its cosm ological 
background (e.g. iBeaslev et ai]|2002l or,more recently, Be kki et al.l 
2008). Differences and similarities between GC and the underlying 
stellar populations have been pointed out in the past. For example. 
Forte, Strom and Strom (1981) found that, in the average, GC ap- 
pear significantly bluer than the stellar halos at all galactocentric 
radii in four Virgo ellipticals. 

In t urn, and after the discovery of GC colour bimodal- 
ity (e.g. IZ epf & Ashman' 1 19931 : lOstrov, Geisler & Fortd Il993h 
iForbes & F orte (2001) noted that the stellar halo colours in the in- 
ner regions of early type galaxies are remarkably similar to those 
of the "red" GC population. 

A more recent comparison between the integrated features of 
the NGC 3923 galaxy halo and its GC population has been pre- 
sented in .Norris et alj(2008) on the basis of Gemini-GMOS spectra 
and direct imaging data. These authors show that the galaxy halo 
has a chemical abundance close to the average value defined by its 
GC and that both systems are coeval. 

A point to be considered in that kind of comparisons, however, 
is that the GC colour statistics are "number weighted" while halo 
colours are, naturally, "luminosity" weighted. Thus, a coincidence 
between GC and stellar halo colours would be expected only if each 
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GC is formed along a given diffuse stellar mass on a constant lumi- 
nosity basis. 

An alternative approach has bee n present ed in 
iForte. Faifer & Geislei|[?200 5) and forte. Faifer & Geisleill(2007ll 
(hereafter FFG07) in a comparative study of the giant ellipticals 
NGC 1399 and NGC 4486 (M87). The main assumption in these 
works is that GC are good tracers of the spatial distribution and 
chemical abundance spread of "diffuse" stellar populations, formed 
along with the clusters, and that the number of GC, TV, per diffuse 
stellar mass unit A/* with chemical abundance Z, is given by 



dN{Z)/dh/L{Z) = 7exp(-5[2-///]) 



(1) 



which is a generalization of the lZepf & AshmM] jl993h T param- 
eter. In those galaxies, the results show that the number of GC per 
diffuse stellar mass increases when metallicity decreases. 

7 and 5 are free parameters that, once the spatial distribution 
of the GC and their chemical abundance scales are determined, al- 
low a fit of the galaxy surface brightness profiles over several tens 
of Kpc. These fits also provide a good approximation to the galacto- 
centric colour gradients as well as of the difference between galaxy 
halo colours and GC mean colours. In turn, the inferred metallic- 
ities for field stars show a remarkable similarity with the observed 
ones in nearby resolved galaxies. 

This paper aims at further testing the modeling described 
above, taking advantage of the largest and most homogeneous GC 
data set avai lable so far in the literature: The Virgo Advanced Cam- 
era Survey jCote etaklliooj and subsequent paper s). In particu- 
lar, w e refer to the galaxy photometry presented by Ferrarese et al.' 
200ft hereafter F2006) and to the GC systems study by Peng et al., 
20061 : hereafter P2006 ). Figure [T] shows the colour difference be- 



tween stellar halo colours and mean GC colours for a galaxy sub- 
sample (see paragraph [5]l from F2006 and P2006. This diagram is 
consistent with the results in Forte et al. (1981) mentioned above. 

A thorough discussion about the characteristics of the GC for- 
mation efficiencies in terms of galaxy brightness, mass, and envi- 
rorunent within th e Virgo cluster has been recently presented by 
IPeng et alj[(2008)| . This work shows that both the number of GC 
per galaxy brightness or galaxy mass units have a non monotonic 
behaviour reaching a minimum for galaxies with stellar masses 
5.1O"'M0 . 

That overall trend, in turn, is comp atible with th e results ob- 
tained by iRhode. Zepf & SantosI ^20051) and later by ISpitler et all 
( l2008l) . in their discuss ion of galaxies more m assive than that mass, 
and also with those bv lMiller & Lotj j2007h in the analysis of GC 
systems associated to less massive dwarf galaxies. 

In this paper we address t he same subject a lthough with some 
differences in comparison with lPeng et alj(2008)i. namely: 

-GC colour histograms are decomposed in terms of "blue" 
and "red" clusters using the FFG07 approach instead of the non- 
parametric fits performed by P2006. 

-Stellar mass to B luminosity ratios are derived on the basis of 
the FFG07 modeling and used to estimate stellar masses within a 
galactocentric radius of ~ 150 arcsecs, the typical coverage of the 
Virgo ACS. 

-Cluster formation is analyzed in the volume space defined by 
GC formation efficiency, galaxy stellar mass, and projected stellar 
mass densities, as well as in terms of chemical abundances. 



2 MODEL ASSUMPTIONS (AND CAVEATS) 

The GC-stellar haloes link given in FFG07 involves several as- 
sumptions, namely: 

a) GC colour bimodality in fact reflects two different ("blue" 
and "red") GC families which have distinct chemical abundance 
and spat ial distributions. Several arguments supp ort this assump- 
tion (see lKundu & Zepj|2007l : IStrader et alj2007l) . 

b) GC colour histograms are usually decomposed via two- 
Gaussian fits that involve five free parameters. Instead, FFG07 
seek for a link betweeen chemical abundance Z and integrated 
colour by means of an empirical colour-metallicity relation (see 
paragraph [3j- Exploring different trial functions leads Xo N{Z) ~ 
exp{Z — Zi) jZs, where Zi is an initial chemical abundance, as the 
simplest approximation for the colour histogram decomposition. In 
this work, as in FFG07, we set Zi — 0.0035.Zq. The histogram 
fit then requires three parameters: the chemical abundance scale 
lenghts Zs{blue), Zs{red) and the relative fraction of clusters be- 
longing to each GC family. 

The case of the blue GC deserves a particular comment in rela- 
tion with the presence of the so called "blue tilt", i.e., these clusters 
seem to become redder with increasing brightness in some galaxies. 
This effect produces a spread of the blue GC colours that FFG07 
identify as a variation of Zi with cluster mass in the case of NGC 
4486. 

We note that the existence of the tilt in NGC 44 8 6, is a subject 
of deb ate in th e literature, e.g. , IStrader et alj (l2006l):lMieske et al.l 
( l2006h , FFG07,'Kundu"(2008?, or lWaters et alJ(2008)L 

c) Blue and red GC families are assumed to be old and formed 
within a relatively narrow period of time. Argu ments supporting 
this assumption can be found, for example, in iBrodie & Straden 
(2006) (see their figure 7). 

Alternatively, some intermediate age globular candidates have 
been repor ted, among the "re d" globulars in other galaxies, 
by Forbes et al.l (12001) (al so see IPierce et"aill2006al . IPierce et"al] 
r2006h : iHempel et al.ll2007l) . 

d) The ratio of the number of GC to diffuse stellar mass with 
the same Z, is assumed to follow the exponential relation given 
in paragraph [T] For typical values of 5, that ratio increases when 
chemical abundance decreases. This is a "phenomenological" as- 
sumption, that allows a good fit to the surface brightness and halo 
colour gradients in NGC 1399 and NGC 4486, and that we confront 
with the Virgo ACS data in this work. 

A caveat in the whole approach concerns to the fact that 
GC statistics refer to clusters that have survived to a number of 
destructive events of dynamical origin l lCharlton & LagunjI 19951: 
iGnedin & O striker 1997), i.e., these GC might not reflect the orig- 
inal cluster populations. For example, both the blue and red GC 
show cores in their projected galactocentric surface densities (al- 
though with widely different spatial scales) that may be the result of 
that depletion process. These cores contrast with the peaked bright- 
ness profiles in the inner region of galaxies. 



3 THE COLOUR METALICITY CALIBRATION 

An em pirical (C— Ti) colo ur versus \Fe/ H]zuj metallicity relation 
(on the lZinn&Wes^l 19841 scale) was presented in FFG07: 

(C - Ti) = 0.94 + 0.068 ([Fe///],„ + 3.5)^ (2) 

This relation was obtained by combining 100 GC in the MW 
and 98 in other three galaxies whose metallicities were derived 
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Figure 1. (g — z) colour difference between galaxy halo colours (from Ferrarese et al. 2006) and the average GC colours in each galaxy (from Peng et al. 
2006). 



2.5 




Figure 2. Integrated (g — z) colour versus chemical abundance from SSP 
models for an age of 12 Gy. Filled circles: .Mar aston (2005)); open circles: 
iBruzual & Chailol (2003). The continuous line is the empirical calibration 
given by FFG07, adopting [Z / H] = [Fe/H]^^ + 0.131, and {g - z) = 
(C - Ti) - 0.20 (see text). 



from Lick indices. This relation is also adopted in this work af- 
ter transforming the [g — z) colours to (C — T i) as discussed in 
par agraph [4| A compar i son w ith SSP models bv lMarastonI ( l2005h 
and iBruzual & Chariot! ( l2003h , for an age of 12 Gy, is shown in 
Figure [2] 



Considering the {g — z) and [B — I) colours in the lMarastonI 

( I2OO5I) models, that leads to: 



(5 



1.00 (B - /) - 0.68 



(3) 



and the empirical (B — J) versus ((7 — Ti ) colours relation for MW 
clusters given in FFG07: 



(C - Ti) = 1.03 {B- I)- 0.43 



(4) 



indicates that (C — Ti) and {g — z) can be linked just by means of 
a zero point difference, a procedure that we describe and adopt in 
the following paragraphs. 



4 PHOTOMETRIC SCALES 

For the sake of homogeneity with FFG07, we first transform the 
g magnitudes given by F2006 to B magnitudes (Johnson system). 
The B luminosities, in turn, will be used to estimate the galaxy 
masses though the [M/L)b ratios derived from fits to the GC 
colour histograms. 

The Virgo ACS survey has a field coverage of 3.4 arcmin, 
which roughly means a galactocentric range of ~ 150 arcsecs. 
Galaxy g magnitudes were then obtained within an elliptical con- 
tour, with a semimajor axis a= 150 arcsecs, and using the Sersic 
best fit parameters given in F2006 (table 3). In turn, integrated 
B magnitudes, within that galactocentric range, were o btained for 
14 gal axies, fainter than B=l 1.0 mag, in common with lCaon et al.l 
( I1990I) yielding: 

~ (5) 



B = 5(150) + 0.28 (±0.19/\/l4) 

a transformation that we adopt in what follows. 

This connection between the 51(1 50) and B magnitudes can 
be also comp a red w ith that derived from the relations given by 
iRodgers et all l l2006l) (based on stellar photometry), and adopt- 
ing a mean colour [g — 2:)=1.3 for the ACS galaxies, leading to 
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Given the GC colour relations in FFG07, this colour transforms 
to (C-ri) = 1.77 and then: 
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Figure 3. (C — Ti) colour histogram for the NGC 4486 globular cluster 
candidates with galactocentric distances from 120 to 360 arcsecs (thin line) 
using the data from FFG07. The thick line corresponds to clusters within 
150 arcsecs from the galaxy center as given by Peng et al. 2006. The binning 
steps are those defined in this last work. Vertical lines show the positions of 
the "blue" and "red" peaks, {g — z) was transformed to (C — Ti ) by adding 
0.29 mags. (see text). 



B = (;(150) + 0.44. The differences between both zero points give 
an idea about the photometric uncertainties involved in this kind of 
analysis. 

The FFG07 modelling delivers (C — Ti) colours that can be 
transformed to — z) just by adopting a colour offset (as discussed 
in the previous paragraph). First we look for the link between the 
(C — Ti ) and (g — z) colour scales from FFG07 and P2006, respec- 
tively, and using the photometry of the well populated GC system 
of NGC 4486 (VCC 1316). Figure[3]shows the (C-Ti) colour his- 
togram for GC within a galactocentric range of 120 to 360 arcsecs 
from FFG07. This diagram also shows the GC histogram given by 
P2006 for the inner regions of the same galaxy, and adopting a shift 



(C - Ti) - (.g - z) (P2006) = 0.29 



(6) 



that agrees with a previous estimate in FFG07. Vertical lines in Fig- 
ure 3 indicate that the position of the "blue" and "red" GC peaks 
are consistent in both statistics (using the same binning steps as 
in P2006) and no variation of these peaks seem detectable with 
galactocentric radius. The sa me feature hold s furthe r out in galac- 
tocentric radius as shown bv' Kundu & Zepj ( l2007h . We note that 
the blue tail of the GC colour distribution is more populated in the 
outer field, suggesting a somewhat bluer lower boundary in colour. 
The significance of this effect is hard to asses since the outer re- 
gion is more likely to be affected by field interlopers (for example, 
unresolved blue galaxies). 

A similar approach was adopted regarding the galaxy inte- 
grated colours. F2006 give (.g - z) = 1.60 for NGC 4486 while 
the integration of the B profile from Caon et al. (1990), combined 
with the (B-R) colour gradient listed by Michard (2000), leads to 
{B — R) = 1.53 within a galactocentric radius of 150 arcsecs. 



(C - Ti) - (.g - z) (F2006) = 0.17 



(7) 



The zero point difference between the F20 06 and P2006 
colour scales is coherent with the results given by ' janz & Liskeil 
12009) in their analysis of the (g-z) colour magnitude relation of 
early type galaxies. 



5 GC COLOUR HISTOGRAM FITS 

P2006 present {g — z) colour histograms for GC systems in 100 
Virgo galaxies (their figure 2). Our analysis is restricted to galax- 
ies with more than ~ 15 GC and at least five colour bins (after 
removing background contamination). Smaller cluster populations, 
or lower definition colour histograms, prevent meaningful fits. 

There are four galaxy fields (P2006 numbers: 43, 61, 79, 93) 
that, being close to galaxies with large CG populations, might be 
likely contaminated by their GC. In turn, galaxy number 39 is the 
only unimodal galaxy with a narrow red peak. This last GC system, 
and those in galaxies 27, 32, 33 and 65, could not be properly fit 
due to noisy histograms or high background contamination. As a 
result of the criteria mentioned above, and after rejecting these last 
objects, our sample finally comprises 63 galaxies listed in Table[T] 

P2006 perform a non-parametric fit to the GC histogram 
which, as mentioned, is different from the approach presented in 
this paper. Here we adopt a Monte Carlo based approach that seeks 
the determination of the chemical abundance scales, Zs{blue) and 
Zs{red), as well as the fraction of clusters in each GC population 
from the GC colour histograms. 

First, a "seed" globular is randomly generated with a given Z 
(controlled by an initial Zs parameter). This value is transformed to 
[Fe/H]zw and to {g—z) colour, by means of our colour-metallicity 
relation. Interstellar colour excesses and gaussian observational er- 
rors (with a dispersion ± 0.05 in {g — z)) are added to each cluster. 
The resulting synthetic histograms are then compared with those 
given by P2006 using the same variable binning steps adopted by 
these authors. 

The Zs values, and the ratio of clusters belonging to the blue 
or red pop ulations, are iterated until a "quality fit" parameter, as 
defined byl Cote et al] l ll998l) . is minimised. 

Throughout this analysis we use an average E(B-V) = 0.028, 
from the colour excesses given by F2006 which, in turn, are based 
on the Schlegel et al. (1998) maps. This seems justified by the un- 
certainties of these maps although, as noted bvl Peng et al.ll(2008ll 
some noise would be expected on the derived blue luminosities aris- 
ing in dust eventualy associated with each galaxy. 

The resulting Zs scales, number of clusters assigned to each 
GC family and the quality of the fit indicator, are given in Table[T] 
Typical uncertainties connected with counting statistics range from 
0.005 to 0.01 for Zs{blue) and 0.05 to 0.15 for Zs{red). 

A sample of histogram fits are shown in Figure |4] for some 
representative galaxies with different numbers of GC (upper row: 
galaxies with more than 250 GC; middle row: with 50 to 100 
GC; lower row: with less than 50 GC). In this figure the fits are 
smoothed histograms adopting a Gaussian colour kernel of 0.03 
mags on the model GC colours. In turn, the histogram fits for the 63 
galaxies listed in Table[T] including the associated statistical count- 
ing uncertainties for each colour bin are displayed in the Appendix. 

Figure |4] shows the same trend noted already by P2006. in the 
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Table 1. Globular clust er colour histogram fit parameters for Virgo ACS galaxies. Col. [1]: Galaxy number from P2006; Col. [2]: VCC galaxy number from 
iBinggeli. Sandage & T ammann 1985; Col. [3]: Number of blue GC; Col. [4]: Chemical abundance scale for the blue GC in Zq uni ts; Col. [5]: Nu mber of 
red GC; Col. [6]: Chemical abundance scale for the red GC in Zq units; Col. [7]: GC colour histogram quality fit parameter (as in'Cote et al.''l998); Col. 
[8]; Fraction of stellar mass associated with the red GC; Col. [9]: Integrated apparent g magnitude (within a=150 arcsecs); Col. [10]: Integrated (g-z) galaxy 
colours from F2006; Col. [11]: Integrated (g-z) galaxy colours from model fit; Col. [12]: Stellar mass to luminosity (B band) ratio from model fit. 
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Table 2. Table 1 - Continued. 
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Figure 4. Sample GC colour histogram fits. Upper row: VCC 1226, VCC1978, VCC 798; middle row VCC 1938, VCC 1303, VCC 1283;lower row: VCC 
1422, VCC 1828, VCC 1407. Upper right numbers are the galaxy identifications in Peng et al. (2006). 



sense that GC colour bi-modality is clearly present in the bright- 
est galaxies and becomes less evident, or is directly absent, in the 
faintest ones. 



6 THE 7 AND 5 PARAMETERS 

This paragraph revises the values of the 7 and 5 parameters given 
in ITGOV on the basis of the following arguments: 

a) Those authors used [Z/H] instead of [FelH\zw (as 
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Table 3. 7 and <5 parameters for NGC 1399, NGC 1427 and NGC 4486. 





NGC 1399 


NGC1427 


NGC 4486 


7 (10~* units) 


0.73 ± 0.05 


0.65 ± 0.03 


0.75 ± 0.03 


5 


1.60 ±0.15 


1.90 ±0.15 


1.75 ± 0.10 



stated in the paper) as the argument of the adopted colour- 
metallicity relation used to obtain integrated model colours. That 
mistake is corrected in this work. As in FFG07, [Z/H] is on the 
iThomas et al. Il l2004i) chemical abundance scale and connected with 
\FelH]^^ through [Z/JJ ] = \F e/H]^u, + 0.131 as derived by 
iMendel. Proctor & ForbesI ( l2007h . 

b) The 5 parameter was determined in FFG07 only trough 
the fit to the observed brightness profiles of NGC 1399 and NGC 
4486. However, that parameter also controls the resulting galacto- 
centric colour gradient. In this paper, we attempted a sim ultaneous 
fit to b oth the brightness and colour gradients, as given in lMichardI 
I2OOO) for these galaxies, aiming at minimizing the orthogonal sum 
of the rms of the brightness and colour fits. 

c) 7 and 5 w ere also derived for NGC 1427 using data from 
iForte et al. I j200lh . This galaxy was observed adopting the (C — 
Ti ) index and using the same reduction and calibration procedures 
reported by FFG07. 

The revised parameters for the three galaxies are listed in Ta- 
ble[3] In particular, the 5 parameters seem comparable (within the 
uncertainties) then suggesting that a similar pattern connects GC 
and stellar haloes in the three systems. We note that the average 
value, S = 1.75 that we adopt hereafter, is higher than those given 
in FFG07 {S = 1.10 and 1.20 ) for NGC 1399 and NGC 4486, re- 
spectively, after correcting the mistake noted in the first item of this 
paragraph. 



7 INTEGRATED COLOURS, STELLAR 

MASS-TO-LIGHT RATIOS AND GALAXY MASSES 



The model assumptions, i.e.: 

dN/dZ = (N/Zs) exp{~Z/Zs) 



(8) 



where N is the number of globular clusters in the blue or red fami- 
lies and Zs their respective abundance scales, and: 



dN/dM4Z) = 7exp(-5 [Z/H] 



(9) 



imply that the amount of diffuse stellar mass with chemical abun- 
dance Z will be given by: 



M,{Z) = 7"' exp(5 [Z/H]) \ exp{-Z/Zs{blue)) + 



L Zs{blue) 

Nred 



exp(-Z/Zs(red))j (10) 



Zs{red) 

Both Nuue and N^ed have distinct behaviours as a function of 
galactocentric radius which can be used to derive surface brightness 
profiles as shown in FFG07. However, in this work, they represent 
the total number of blue or red GC within the ACS fields. 

The B luminosity associated to Mt{Z) was estimated by 
adopting a reference mass-to-luminosity ratio: 



{M/L)b = 3.71 + {[Z/H] + 2.0)' 



(11) 



with {M/L)b = 3.71 for [Z/H] < -2 and {M/L)b = 13.6 
for [Z/H] > 0.5, which gives a good approximation (within ±7 



percent) to the SSP models bv lWorthevlll994h for an age of 12 Gy 
and a Salpeter luminosity function. A comparison with other mod- 
els gives an idea about the uncertainty of this ratio. For example, 
models by Maraston (2005), for the same age and luminosity func- 
tion, shows an overall agreement within 10 percent with Worthey's 
except at the lowest abundance where the ratio is about 24 percent 
larger. 

The total stellar mass, integrated B luminosity and compos- 
ite [M/L)b ratio can be found by integrating M,{Z) dZ and 
Mt{Z) {AI/L)g^ dZ, respectively, where we adopted Zi = 
0.0035 Zq and Zu ~ 4.0 Zq as lower and upper integration limits. 

In what follows, we tentatively set S = 1.75 as representative 
for all galaxies (see paragraph [SJ. With this assumption, the inte- 
grated mass-to-light ratios, as well as the total mass fraction in a 
given stellar population, will only depend on the Nuue, Zs{blue), 
Nred and Zs{red) values obtained from the GC colour histogram 
fits and are independent of the 7 parameter. This also holds for the 
integrated galaxy colours that, for each M{Z), depend only on the 
chemical abundance through the adopted colour-metallicity rela- 
tion and on the value of 5, which determines the amount of stellar 
mass with a given abundance. 

Table[T]lists the stellar mass fraction associated to the red GC 
population, F, the integrated g(150) magnitude, the {g — z) galaxy 
colours from F2006, the {g — 2) colours derived from the model 
and, finally, the inferred stellar mass to B luminosity ratio. 

In the case of M87, the adoption of the parameters given in 
Table[T]leads to (g — z) = 1.55 , which is 0.05 bluer than the value 
quoted by F2006. This difference can be removed by increasing the 
Zs{red) scale by 1.25, i.e., a fraction of the stellar mass associated 
with the red GC, in the inner 150 arcsecs, seem to have formed 
with Z values somewhat higher than in these clusters. We note that 
this is not the case at galactocentric radii larger than 120 arcsecs 
where the integrated colours of both NGC 1399 and NGC 4486 
can be matched with the Zs(red) derived from the histogram fits 

(as shown in FFG07). 

An argument in the same sense can be found in lBeaslev et al.l 

( l2008l) who show that, in the inner regions of the resolved galaxy 
NGC 5128, field stars have upper [Z/H] values some 0.2 dex 
higher than those of the red GC. On the basis of these considera- 
tions we adopted a chemical abundance spread equal 1.25 Zs(red) 
in the estimate of the integrated colours of all our sample galaxies. 

The residual colours (observed minus model) as a function of 
the integrated (g — z) colours given by F2006 are shown in Figure|5] 
(left panel). The expected impact on the calculated (g — z) colours, 
arising the uncertainties of the chemical abundance scales referred 
to in paragraph [5] ranges from 0.02 to 0.06 mags. The larger dis- 
persion in this diagram then probably arises in observational errors 
that become more important for the bluer (and fainter) galaxies. 

This diagram shows that the model is able to approximate 
the observed galaxy colours although an average {g — 2) residual 
of 0.07 mags remains. This residual, as shown in Figure |5] (right 
panel), correlates marginally with the distance to the centre of the 
Virgo cluster (that we tentatively identify with NGC 4486) in the 
sense that, inner galaxies, appear slightly redder than the colours 
predicted by the model. We note that galaxies number 73 (VCC 
21; not included in Figure 5a) and number 84 (VCC 1440), exhibit 
(g — 2) colours that are significantly bluer than those inferred from 
their GC and possibly denoting the presence of younger stellar pop- 
ulations. 

Stellar mass-to-B luminosity ratios as a function of stellar 
mass are depicted in Figure |6] (left panel). The v ertical line in this 
diagram shows the range in {M/L)b obtained bv lNapolitano et all 
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Figure 5. Left panel: (g — z) colour differences between the observed halo colours (from Ferrarese et al. 2006) and those obtained from modelling as described 
in this paper. Right panel: Stellar haloes (g — z) colour residual (observed minus model) as a function of the angular distance to M87. The straight line suggests 
a marginal trend of colour residuals with position within the Virgo cluster. 



( l2005h . These authors present a comparison between ratios derived 
through dynamical and SSP models that, according to the last dia- 
gram, are in very good agreement with the values obtained in this 
paper. We stress that the relation between the {M/L)b values and 
the blue absolute magnitude of the galaxies (listed in Table 1 and 3, 
respectively) show good agreement w ith the analysis of the SDSS 
data given by iKauffmann et alj Jioos l). In particular, those values 
define an upper envelope to their relation suggesting that the VCC 
galaxies in our sample are in fact the oldest for a given galaxy 
brightness. 

Total masses within a=150 arcsecs were derived by transform- 
ing g{150) magnitudes to B absolute magnitudes (as explained in 
paragraph Ht, adopting adistance modulus to Virgo, {V — Mv)o = 
31.1 from lTonrv et^feOOlh . and the [M /L) b ratios given in Ta- 
ble|4](column 6). This table also lists the absolute magnitudes Afg 
and angular distance of each galaxy to NGC 4486 (in degrees). 

The d erived stellar masses M, are compared with those given 
bv lPeng e t al. (2008) in Figure[6](right panel). This comparison ex- 
cludes the eight brightest galaxies since, being extended objects, 
they have a considerable fraction of their stellar masses beyond 
a=150 arcsecs. The straight line 



Zo(/(M,)peng et al. = 1-0 log{A'Q ~ 0.25 



(12) 



indicates that our masses are a factor about 1.8 larger This differ- 
ence arises since our reference mass to luminosity ratio (eqn. 11) 
adopts a Salpeter stellar mass distribution while Peng et al. (2008) 
models stand on that derived by IChabrieili(2003)l . The GC forma- 
tion efficiency discussed later, however, will keep the same depen- 
dence on stellar mass as the slope of the relation between both mass 
estimates is not significantly different from unity. 



8 CHEMICAL ABUNDANCES 

The Zs{blue) and Zs{red) parameters (from Table[TJ are depicted 
as a function of total stellar mass in Figure [7] (left panel). This di- 
agram shows that the chemical abundance scale of the blue GC 
keeps mostly below 0.05 Zq. Five galaxies with Af, smaller than 



10" Mq (numbers 50, 91, 62, 67, 95), among the 18 galaxies 
that exhibit unimodal histograms, show rather high Zs(blue) val- 
ues. Their position in that diagram suggests they might be in fact 
bimodal, i.e, include a red GC component not much different in 
chemical abundance scale from that of the blue GC and then broad- 
ening the histogram. 

For stellar masses larger than 10^" Mq there seems to be a 
mild increasing trend of Zs{blue) with galaxy mass as shown in 
Figure|7] right panel. This diagram also shows a linear least squares 
fit: log{Zs{hlue)) = -0.07 + 9.10"^ log{M,). This fit does not 
include the five galaxies mentioned before as is restricted to galax- 
ies with stellar masses larger than 10^'' Mq. 

This trend is coherent with IStrader. Brodie & ForbesI i2004h 
(see their figure 1) who find an increase of the mean [V—I) colours 
of blue GC with increasing galaxy brightness. 

Figure |7] also shows that the GC Zs{red) parameters (filled 
dots) clearly correlate with galaxy mass. This result supports the 
idea that red GC can reach higher metallicities in the deeper poten- 
tial well of the more massive systems. The relatively large spread 
of this relation, in turn, may be reflecting different degrees of het- 
erogeneity connected with the formation of the red GC and their 
associated diffuse stellar population. 

Blue luminosity weighted abundances [Z/H], derived from 
the model are shown as a function of stellar mass and abso- 
lute B magnitude in Figure [S] and listed in Table|4] (column 5). 
The first of these figures also includes a straight line with a 
slope d\Z/H]/dlog( M,) = —0.20 , comparable to values found 
in the literature (se e iLa Barbera et alj |2008| , iKodama et al. 
iThomas et ai]|2005l) . and gives a good approximation covering the 
upper two dex in stellar mass. However, both diagrams show a 
change in the slopes at a stellar mass A/, ~ 10^" A/© or M{B) ~ 
— 17.5. We note below that stellar mass, the GC colour distribu- 
tion becomes predominantly unimodal, i.e., the red GC compo- 
nent seems very small or is directly absent. The referred mass then 
seems then a key parameter regarding the type of stellar population 
present in a given galaxy. 
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Table 4. Chemical abundances, structural paramteters and globular cluster formation efficiency parameters for Virgo ACS galaxies. 
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Figure 6. Left panehStellar mass to Blue luminosity ratio as a function of total stellar mass derived in this paper for 63 galaxies in the ACS Virgo Survey. 
Right panel: Comparison between the stellar masses derived in this paper with those given by Peng et al. (2006). The straight line (with unity slope) is the 
approximation given in the text. 
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Figure 7. Left panel: Chemical abundance scales Zs(blue) and Zsired) for the blue and red GC (open and filled dots, respectively) as a function of galaxy 
stellar mass. Blue GC show a mild increase with mass consistent with previous results (e.g. Strader & Brodie 2007). Right panel: Zs{blue) values, for galaxies 
more massive thanM* IQI^Mq, as a function of galaxy stellar mass on an expanded scale. The straight hne is a least square fit to the data points. 



9 THE GC FORMATION EFFICIENCY-GALAXY 
MASS-PROJECTED MASS DENSITY SPACE 

The GC formation efficiency t can be defined as a function of the 
total numer of GC and total stellar mass A/, or, discriminating both 
GC famihes, as tred = Nred/M* and ttiue = Nbiue/Mt. The be- 
havior of these three parameters as a function of the galaxy struc- 
tural parameters is analyzed in what follows. 

Figure|9](left panel) shows the (logarithmic) volume space de- 
fined by stellar mass, projected stellar mass density a (within the 
effective radius re and in solar masses per sq. Kpc), and total GC 
efficiency, t (all values listed in Table|4j. These parameters were 
derived by adopting our stellar mass estimates, the galaxy effective 
radii from F2006, and the number of GC listed in Table[T] 

This figure shows that most galaxies with stellar masses be- 



low ~ lO^'^ Mq define a thick plane. In turn, eight of the brightest 
galaxies appear clearly detached and above that plane. These last 
galaxies exhibit intermediate to low projected mass densities and, 
also, low central surface brightness co mpared to the Sersic p rofile 
that fits the outer galaxy profile (see 'Kormen dv et ai]|2008 ). We 
point out that the t parameters listed in Table|4]for these galaxies 
are only approximate values. This is due to the fact that both GC 
families have cores in their spatial distribution (i.e. the projected 
GC density profiles are shallower than the galaxy brigtness profiles) 
leading to low t ratios in their inner regions. On the other side, as 
noted before, their total masses include a considerable fraction be- 
yond a=150 arcsecs. A tentative correction for this last effect would 
increase their tabulated log{a) values by ~ 0.25. 

A least square fit to the t values (in 10~^ units) and assuming 
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Figure 8. Left panel: Blue luminosity weighted chemical abundance as a function of total stellar mass. The straight line as a slope of -0.2, comparable to 



values reported in the literature. A significant change in that slope appears for masses smaller than Mt 
(within a= 150 arcsecs) versus chemical abundance (B luminosity weighted). 
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Figure 9. Left panel: Globular cluster efficiency (total number of clusters per stellar mass unit) as a function of stellar mass and projected stellar mass density 
(a). Filled dots lie below the plane. Dots with a cross belong to the most massive galaxies. Right panel: Edge on view of the plane depicted in the previous 
figure, for galaxies with total stellar mass smaller than 10^^ ^Q- 



no errors on effective radii or projected mass density, yields, for 49 
galaxies with stellar masses smaller than 10^^ Mq: 

logit) = -0.70(±0.05) log{M*) + 0.11(±0.04) log{a) + 

6.75(±0.04) (13) 

with a fit rms of ±0.17 in log(t). The fit residuals in log(t) show 
no correlation with galaxy mass, integrated colour, morphology or 
position within the Virgo cluster. 

An "edge-on" view of the plane is depicted in Figure[9](right 
panel) excluding the eight brightest galaxies in Table[T] 

The dependency with projected mass density becomes more 
evident if only the red GCs (in 31 bimodal galaxies), are included 
in the fit: 



log{Ued) = -0.65(±0.16) /off(M.) + 0.21(±0.09) log{a) + 

4.99(±0.05) (14) 

with a fit rms of ±0.21 in log{t,.ed)- 

On the other side, blue GC in 49 galaxies (3 1 with bimodal and 
18 with unimodal GC systems) show a steeper slope with stellar 
mass but no detectable dependence with projected mass density: 

logihiue) = -0.91(±0.07) Zog(Af.) + 0.02(±0.05) log{a) + 

9.37(±0.04) (15) 

with a fit rms of ± 0.19 in log{tbiue)- 

This lack of dependency may be expected as the dominant 
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stellar mass in bimodal GC galaxies appears associated with the 
red GC population according to Table[T] 



10 CONCLUSIONS 

The analysis presented in this work suggests that GC in early type 
galaxies are good tracers of their stellar populations. We stress that, 
the proposed GC formation efficiency dependence with chemical 
abundance, leads to a good approximation of the galaxy integrated 
colours and is consistent with a constant 5 parameter suggesting a 
common pattern for the GC-field stars connection. 

We also note, in particular, that this approach is able to explain 
the stellar halo-mean GC colour differences without implying def- 
init ely distinct format ion histories for field stars and GC, although 
see IPeng et alj(2008il for a different interpretation. 

In turn, the inferred galaxy (M / L) b ratios compare well with 
the range of values found in the literature and seem consistent with 
old ages and small formation time spreads. 

The on-set of GC colour bimodality seems to occur at masses 
larger than A/, ~ 1O^"A/0, i.e., galaxies with lower masses only 
give rise to blue GC. Interestingly, Dekel & Woo ( 2003 ) suggest 
that galaxies with masses smaller than « 'i.lQi^^ Mq are efficiently 
cleaned of interstellar matter as a result of stellar and SNe winds 
leading to an interruption of the star formation process. 

The two dominant GC families have rather distinct character- 
istics. On one side, "blue" globulars can be characterized by a very 
low chemical abundance spread suggesting that some kind of major 
common event may have played a role in their "synchronous" for- 
mation and relatively homogeneous chemical composition. The re- 
ionization of the Universe has b een put forwar d as possible respon- 
sible for this characteristic (e.g. lSantosll2003[) although the impact 
of re-ionization on galaxy and galaxy-cluster scales should also de- 
serve further analysis. In any case, the increase of Zs(hlue) with 
galaxy mass lends support to the idea th at blue GC in fact "know" 
about the galaxy they are related to (e.g. I S trader. Brodie & Forbe3 
l2004h . 

The very narrow metallicity scale of the blue GC is compatible 
with their formation in low effective-yield environments character- 
ized by important mass loss leading to a rapid supression of field 
star formation (e.g. [Harris et"al]|2006h and to high GC formation 
efficiencies. 

In turn, red clusters exhibit an abundance spread that becomes 
larger for more massive galaxies and appear associated, when GC 
bimodality is observable, with the dominant stellar mass compo- 
nent. This is coherent with the idea that more massive systems ex- 
hibit higher mean chemical abundances. The large dispersion of 
the relation between the Zs{red) scale and galaxy mass may be 
also indicating some diversity in the formation history of this com- 
ponent (e.g. different numbers of mergers and their effects on gas 
removal). 

The overall increase of chemical abundance with mass derived 
in this work, is consi stent w ith the GC formation modelling pre- 
sented bv lBekki et al.l ( l2008l) . These models imply a formation time 
spread of roughly 3 Gy (i.e. ±1.5 Gy around a mean value of about 

11 Gyr. 

We also stress that the broad behaviour derived for t with 
galaxy mass is coherent with the starting hypothesis in this work, 
i.e., that the number of clusters per diffuse stellar mass unit in- 
creases with decreasing metallicity (and hence with decreasing stel- 
lar mass according to the stellar mass-metallicity relation). 

In more detail, galaxies with stellar masses lower that M, ~ 



10^^ Mq define a thick plane in the logarithmic stellar mass-tj-t 
space where the GC formation efficiency seem dependent on both 
galaxy mass and projected surface mass density. Such a dependence 
becomes more evident when only r ed GC are taken into con sidera- 
tion. This result resembles those bv lLarsen & Richtlerl ( 12000 ) who 
suggest a dependence of the formation of young and massive stellar 
cluster with projected mass density. The lack of such a dependence 
in the case of the blue GC reinforces the idea that these clusters 
have genetic differences with the red ones. 

Remarkably, the most massive galaxies in the sample appear 
clearly detached and above that plane. Naively, this is the expected 
situation for galaxies suffering mergers, although a number of dy- 
namical effects that may have an impact on the t parameter still 
remain to be considered before adopting this conclusion. For exam- 
ple, in the case of (dry) merging two galaxies with similar masses 
and GC populations, the composite t parameter would remain un- 
changed while the total mass would be doubled. These massive 
galaxies also exhibit large effective radii and intermediate pro- 
jecte d stellar mass densiti es, a signature of major merger events 
(e.g. lShin & Kawatj|2009l) . 

A tentative conclusion, that will deserve further analysis, sug- 
gests that galaxies with stellar masses below that value, and located 
on the "thick" plane described in this paper, may have suffered 
a relatively small number of disippationless merging events along 
their life. This conclusion also g ets support from CDM based mod- 
elling (e.g. lDe Lucia et alj2007h showing that less massive galaxies 
may have had a low number of galaxy-progenitors in their forma- 
tion history. 



ACKNOWLEDGMENTS 

This work was supported by grants from La Plata National Univer- 
sity, Agenda Nacional de Promocion Cientffica y Tecnologica, and 
CONICET, Argentina. JCF acknowledges Dr. Michael West hospi- 
tality at ESQ during the last stages of this paper. Useful comments 
and advise from an anonymous Referee are also acknowledged. 



REFERENCES 

Beasley M. A., Baugh CM., Forbes D.A., Sharpies R.M., Frenk 

C. S., 2002, MNRAS, 333, 383 
Beasley M. A., Bridges T, Peng E., Harris W.E., Harris G.L.H., 

Forbes D.A., Mackie G, 2008, MNRAS, 386, 1443 
Bekki, B. Yahagi H., Nagashima M., Forbes D.A., 2008, MNRAS 

387, 1131 

Binggeh, B., Sandage A., Tammann G.A., 1985, AJ, 90, 1681 

Brodie J.P, Strader J. 2006, AR&A 44, 193 

Bruzual G., Chariot S., 2003, MNRAS 344, 1000 

Caon N., Capaccioh M., Rampazzo R., 1990, MNRAS, 242, 24 

Chabrier G., 2003, PASP 115, 763 

Charlton J.C., Laguna P, 1995, ApJ, 444, 193 

Cote P Marzke R. O., West M. J., 1998, ApJ, 501, 554 

Cote P., Blakeslee J.P, West M. J., Ferrarese, L., Jordan A., Mei 

S., Merritt D., Milosavljevic M., Peng E.W., Tonry J. L., West 

M. J., 2004, ApJS, 153, 223 
De Luci'a G., Springel V., White S.D.M., Croton D., Kauffman G., 

2006, MNRAS, 366, 499 
Dekel A., Woo J., 2003, MNRAS, 368, 2 
Eggen O., Lynden Bell D., Sandage A.R., 1962, ApJ, 136, 748 



The Globular clusters-stellar haloes connection in early type galaxies 13 



Ferrarese L., Cote P.Jordan A., Peng E.W., Blakeslee J.P, Piatek 

S., Mei S., Merritt D., Milosavljevic M., Tonry, J. L., West M. J., 

2006, ApJS, 164, 334 
Forbes D. A., Beasley M.A., Brodie J.P, Kissler-Patig M., 2001, 

ApJ, 563, 143 
Forbes D. A., Forte, J.C., 2001, MNRAS, 322, 257 
Forte J.C., Strom S. E., Strom K., 1981, ApJ, 245, L9 
Forte J.C., Geisler D., Ostrov PG., Piatti A.E., Gieren W., 2001, 

AJ, 121, 1992 

Forte J.C., Faifer R, Geisler D., 2005, MNRAS, 357, 56 

Forte J.C., Faifer F, Geisler D., 2007, MNRAS, 382, 1947 

Gnedin O.Y., Ostriker J.P, 1997, 474, 223 

Harris W.E., van den Bergh S., 1981, AJ, 86, 1627 

Harris W.E., Whitmore B.C., Karakla D., Okon W., Baum W.A., 

Hanes D.A., Kavelaars D.J., 2006, ApJ, 636,90 
Hempel M., Kissler-Patig M., Puzzia T.H., Hilker M., 2007, 

A&A, 463, 493 
Janz J., Lisker T, arXiv 0903.5301vl 
Kauffmann, G. et al. 2003, MNRAS 341, 33 
Kodama T. et al. 1998, A&A, 334, 99 

Kormendy J., Fisher D.B., Cornell M.E., Bender R., 2008, arXiv 
0810.1681 

Kundu A., Zepf, S.E., 2007, ApJ, 660, L109 
Kundu A., 2008, AJ 136, 1013 

La Barbera F, Busarello G., Merluzzi P., De la Rosa I., Coppola 

G., Haines CP, 2008, ApJ (accepted) astro-ph 0807.3829 
Larsen S. S., Richtler T, 2000, A&A, 354, 836 
Maraston C, 2005, MNRAS, 362, 799 

Mendel J. T, Proctor R. N., Forbes D. A., 2007, MNRAS, 379, 
1618 

Michard, R., 2000, A&A 360, 85 

Mieske, S. et al.,2006, ApJ, 653,193 

Miller B.W., Lotz J., 2007, ApJ, 670, 1074 

Napolitano N.R., Capaccioli A.J., Romanowsky A.J., Douglas 

N.G., Merrifield M.R., Kuijken K., Arnaboldi M., Gerhard O., 

Freeman K.C., 2005, MNRAS, 357, 691 
Norris M.A., Sharpies R. M., Bridges T., Gebhardt K., Forbes D. 

A., Proctor R., Faifer F, Forte J.C., Beasley M.A., Zepf S. E., 

Hanes D.A., 2008, MNRAS 385, 40 
Ostrov PG., Geisler, D., Forte, J.C., 1993, AJ, 105, 1762 
Peng E. W., et al., 2006, ApJ, 639, 838 
Peng E. W., et al., 2008, ApJ, 681, 197 

Pierce M., Beasley M. A., Forbes D. A., Bridges T., Gebhardt K., 

Faifer F. R., Forte J. C, Zepf, S. E. Sharpies R., Hanes D. A., 

Proctor R., 2006a, MNRAS 366, 1253 
Pierce M., Bridges T., Forbes D. A., Proctor R., Beasley M. A., 

Gebhardt K., Faifer F R., Forte J. C, Zepf S. E., Sharpies R., 

Hanes D. A., 2006b, MNRAS 368, 325 
Rhode K.L., Zepf S.E., Santos M.R., 2005, ApJ 630, L2 
Rodgers C. T, Canterna R., Smith J.A., Pierce M.J., Tucker D.L., 

2006, A J, 132, 989 
Santos M.R., 2003, in Extragalactic Globular Cluster Systems, ed. 

M. Kissler-Patig (New York: Springer Verlag) 
Schlegel D., Finkbeiner D., Davis M., 1998, ApJ, 500,525 
Searle L., Zinn R., 1978, ApJ, 205, 357 
Shin M-S.,KawataD., 2009, ApJ, 691, 83 

Spitler L.R., Forbes D.A., Strader J., Brodie J.P, Gallagher J.S., 

2008, MNRAS, 385,361 
Strader J., Beasley M.A., Brodie, J.P, 2007, AJ, 133, 2015 
Strader J., Brodie J. P, Forbes D. A., 2004, AJ, 127, 295 
Strader J., Brodie J. P, Forbes D. A., 2004, AJ, 127, 3431 



Strader J., Brodie J. P, Spitler M., Beasley M.A., 2006, AJ, 132, 
2333 

Thomas D., Maraston C, Korn A., 2004, MNRAS, 351, L19 
Thomas D., et al., 2005, 2005, ApJ, 621, 673 
Tonry T, Dressier A., Blakeslee J.P, Ajhar E.A., Fletcher A.B, 
Luppino G.A., Metzger M.R., Moore C.B., 2001, ApJ, 546, 681 
Waters C.Z., Zepf S.E., Lauer T, Baltz E.A., arXivOSl 1.0391 
Worthey G., 1994, ApJS, 95, 107 
Zepf, S.E., MNRAS 264, 611 
Zinn R., West M. J., 1984, ApJS, 55, 45 
Zepf S. E., Ashman K. M., 1993, MNRAS, 264, 611 



APPENDIX:GC COLOUR HISTOGRAMS 

GC colour histogram fits for 63 galaxies (*) in the Virgo ACS. 
Heavy black lines are the (background corrected) observed his- 
tograms. Red lines show the histograms that give the best fits in 
terms of the quality indicator (see Table[T). Dots show the counting 
uncertainties for each colour bin. Note that Figure 4, in this paper, 
only shows nine representative histogram fits. 
(*) In the print version of the journal, the figures for only the first 12 
galaxies are included. The figures for all 63 galaxies can be found 
in the online version. 
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